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I.  BRIEF  OUTLINE  OF  RESEARCH  FINDINGS 


During  the  last  year  we  have  worked  on  problems  of  wave  propagation  in  random 
media,  reflection  of  waves  from  rough  surfaces,  effective  parameters  of  heterogeneous  ma¬ 
terials,  inverse  scattering  theory,  acoustoelasticity,  Hill’s  equation  and  some  other  topics. 
We  shall  describe  some  of  this  work  in  this  section.  Most  of  the  rest  of  our  completed  work, 
which  is  not  described  here,  is  contained  in  the  research  papers  which  are  listed  in  Section 
II.  Some  of  the  articles  listed  there  have  been  published  already,  others  have  been  accepted 
for  publication  and  arc  now  in  press  and  some  have  been  submitted  for  publication  but 
have  not  yet  been  accepted.  The  status  of  each  paper  is  indicated  after  its  title,  and  the 
references  to  the  published  papers  arc  also  given. 

Concerning  waves  in  random  media.  Dr.  Graham  Eatwcll  and  Professor  J.B.  Keller 
have  treated  w.avcs  in  a  medium  of  randomly  distributed  discrete  scatterers.  They  have 
used  a  method  like  the  path  integral  procedure  used  by  Dashen  to  Diul  the  moments  of  a 
Geld  in  a  continuous  random  medium.  A  paper  on  this  work  is  in  preparation. 

Dr.  .lohn  G.  Watson  and  Professor  Keller  h.ivc  published  two  papers.  Nos.  83  and 
112,  on  reflection  of  waves  from  rough  boundaries.  They  arc  now  completing  a  companion 
p.aper  on  rough  intcrfiaccs.  The  novel  features  arc  the  calculation  of  second  moments  of 
the  sc.attcrcd  Geld,  and  the  relation  of  the  results  to  Twersky's  method. 

The  effective  viscosity  of  a  periodic  suspension  .and  the  effective  elasticity  tensor  of 
a  periodic  composite  have  been  calcuhatcd  by  Dr.  Kevin  C.  Nunan  and  Professor  Keller 


iu  papers  Nos.  100  and  101.  The  new  features  are  that  all  components  of  the  viscosity 
and  elasticity  tensors  are  found  at  all  concentrations  up  to  close  packing.  Related  work  on 
effective  conductivity  b  being  done  by  Mr.  John  Nevard  and  Professor  Keller.  So  far  the 
reciprocal  theorem  has  been  extended  to  anisotropic  media. 

In  inverse  scattering  theory  Dr.  Margaret  Cheney  has  developed  the  two  and  three 
dimensional  theory  in  a  number  of  papers  (86,  99,  125),  and  applied  it  in  several  others 
(104,  123,  124)  to  problems  which  are  more  realistic  than  many  which  were  treated  pre¬ 
viously.  Dr.  John  Fawcett  has  also  developed  some  theoretical  results  on  stability  (116) 
and  on  inversion  of  spherical  means  (121)  as  well  as  several  more  applied  results  (117, 
118,  120).  Professor  Stephanos  Venakidcs  has  employed  the  inverse  scattering  method  to 
solve  the  Korteweg-de  Vries  equaion  in  the  zero  dispersion  limit  (in  paper  No.  106)  and  in 
a  manuscript  in  preparation.  Mr.  William  Boyse  and  Professor  Keller  are  analyzing  the 
inverse  problem  for  clastic  media  by  the  Born  approximation. 

Dr.  Luis  L.  Bonilla  has  studied  the  effective  elastic  constants  of  aggregates  of  crys¬ 
tallites,  each  with  cubic  symmetry,  by  using  invariance  and  symmetry  methods  (115).  He 
and  Professor  Keller  have  also  treated  wave  propagation  and  the  acoustoclastic  elfect  in 
such  media  (119). 

Dr.  Michael  I.  Weinstein  and  Professor  Keller  have  successfully  analyzed  Hill’s  equa¬ 
tion  with  a  large  potential  and  found  the  behavior  of  the  stability  bands  in  paper  122  and 
in  a  paper  in  preparation. 


III.  ABSTRACTS  OF  MANUSCRIPTS  SUBMITTED  DURING  REPORT  PERIOD 


Newton’s  Second  Law,  by  J.B.  Keller 

Newton’s  second  law  of  motion  is  usually  written  in  the  form  /  =  ma.  Here  o  denotes 
the  acceleration  of  a  point  particle,  m  denotes  the  nmss  of  the  particle,  and  f  denotes  the 
force  acting  on  the  particle.  We  wish  to  elucidate  the  empirical  content  of  this  law,  and  in 
particular  to  indicate  whether  it  emhodics  a  definition  of  force,  a  definition  of  mass,  or  a 
definition  of  both  force  and  mass. 

Effective  Elastic  Constants  of  Polycrystallinc  Aggregates,  L.L.  Bonilla 

A  method  is  presented  for  the  determination  of  the  effective  elastic  constants  of  a  trans¬ 
versely  isotropic  aggregate  of  crystallites  with  cubic  symmetry.  The  results  obtained  gen¬ 
eralize  those  given  in  the  literature  for  the  sccon<l  and  third  order  clastic  constants.  In 
addition,  the  second  iiionicnts  and  the  binary  angular  correlations  of  the  second  order 
stiffnesses  are  olitained.  It  is  also  explained  how  these  moments  can  be  used  to  find  the 
two-point  correlations  of  the  clastic  constants. 

On  the  Stability  of  Inverse  Scattering  Problems,  J.A.  Fawcett 

Some  aspects  of  the  time-dependent  inverse  scattering  i>roblcm  arc  discussed.  Then  some 
simple  stability  results  arc  presented.  Several  analytical  and  nniiicrical  cxam]>les  arc  given 


to  illustrate  them. 


Two  Dimensional  Modelling  and  Inversion  of  the  Acoustic  Wave  Equation  in 
Inhomogeneous  Media,  J.A.  Fawcett 

In  this  paper  we  outline  a  fast  and  eflicient  method  for  the  finite  difference  calculation  of 
solutions  to  the  two  dimensional  acoustic  wave  equation  in  an  inhomogeneous  medium. 
We  then  implement  our  forward  modeling  scheme  in  an  inversion  procedure  to  determine 
the  two  dimensional  velocity  variation.  A  numerical  example  is  given. 

Three  Dimensional  Ray  TVacing  and  Geophysical  Inversion  in  Layered  Media, 
J.A.  F.awcctt  and  II. B.  Keller 

In  this  paper  the  problem  of  finding  seismic  rays  in  a  three  dimensional  layered  medium 
is  examined.  The  “layers”  are  separated  by  arbitrary  smooth  interfaces  that  ciin  vary  in 
three  dimensions.  The  endpoints  of  each  ray  and  the  seqticncc  of  interfaces  it  encounters 
are  specified.  The  problem  is  formulated  as  non-linear  system  of  equations  and  efficient, 
accurate  methods  of  solution  arc  discussed.  An  importtant  application  of  ray  tracing  meth¬ 
ods,  which  is  discussed,  is  the  non-linear  least  squares  cslimation  of  medium  ]>nramcters 
from  observed  travel  times.  In  addition  the  “type”  of  each  ray  is  talso  delermincd  by  the 
least  squares  process  —  this  is  in  effect  a  deconvolution  j)rocerliirc  similar  to  that  desired 
in  seismic  exploration.  It  enables  more  of  the  measured  data  to  be  used  without  filtering 
out  the  miillijile  redections  that  are  not  pure  P-waves. 

Acousloc'lnstic  Effects  and  Sound  Wave  Propagation  in  llctcrogeous  Anisotropic 
Materials,  L.L.  Bonilla  and  J.B.  Keller 

A  theoretical  analysis  of  acoustoclastic  clfccts  is  prc.senlcd.  It  is  based  upon  the  theory  of 
sound  wave  propagation  in  a  stressed  heterogeneous  anisotropic  elastic  medium  composed 


of  grains.  The  ciTcct  of  residual  stress  is  included,  and  shown  to  be  different  from  that 
of  applied  stress.  The  statistics  of  grain  orientation  and  grain  correlation  are  taken  into 
account.  The  acoustoelastic  coefficients  and  the  effects  of  dispersion,  attenuation  and 
texture  of  the  medium  arc  determined. 

Tomographic  Reconstruction  of  Velocity  Anomalies,  J.A.  Fawcett  and  R.W.  Clay¬ 
ton 

An  approximate  inversion  formula  is  proposed  for  the  reconstruction  of  slowness  anomalies 
in  a  known  tlcptli  varying  background  field.  The  data  arc  observed  travcitime  perturbations 
for  reflections  from  a  known  plan<*ir  reflector.  The  limitations  of  the  formula  arc  discussed 
and  numerical  examples  arc  given. 

Inversion  of  N-dimcnsional  Spherical  Averages,  J.A.  Fawcett 

In  this  paper  we  give  a  simple  inversion  formula  for  determining  the  even  part  of  a  function 
from  its  averages  over  the  surfaces  of  spheres  in  IR”,  conlred  on  a  hyperplanc  y  =  0.  We  first 
give  the  formula  in  Fourier  space  and  then  in  physical  space  its  a  generarued  convolutional 
back- projection. 


Hill's  Equation  with  a  Large  Potential,  M.l.  Weinstein  and  J.ll.  Keller 

We  obtain  asymptotic  e\|)aii.sions,  for  large  values  of  the  parameter  A,  of  the  stability 
boundaries,  the  stability  band  widths,  the  Floquet  multipliers  and  the  solutions  of  Hill’s 


equation 


The  potential  ^(x)  is  assumed  to  be  periodic  and  to  have  a  unique  global  minimum  within 
each  period,  at  which  q"  >  0.  The  results  for  the  stability  band  widths  show  that  they 
decay  exponentially  with  A  as  A  increases.  These  results  generalize  those  for  symmetric 
potentials  due  to  Harrell,  and  that  for  the  Mathieu  equation  due  to  Mcixncr  and  Schafke. 

Inversion  of  the  2.5-d  Acoustic  Equation,  M.  Cheney,  S.  Coen  and  A.  Weglein 
An  inverse  acoustic  scattering  theory  and  algorithm  is  presented  for  the  reconstruction 
of  a  two  dimensional  inhomogeneous  acoustic  medium  from  surface  measurements.  The 
measurements  of  the  surface  pressure  due  to  a  luarmonically  oscillating  surface  point  source 
at  two  arlutrary  frequencies  allows  the  separate  reconstruction  of  the  density  and  velocity 
of  the  subsurface.  The  solution  of  the  inverse  acoustic  scattering  problem  is  obtained  by  a 
transformation  to  an  inverse  quantum  scattering  problem  in  two  dimension,  the  solution 
of  which  has  been  obtained  recently  by  Cheney,  who  was  inspired  by  Newton’s  solution  in 
three  dimensions.  The  solution  presented  is  a  generalization  of  Coen’s  and  Stickler  and 
Ueift’s  solutions  of  the  three  dimensional  acoustic  problem  for  a  plane  stratified  acoustic 
medium. 

The  Connection  Between  Time-  and  Frequency-domain  Three  Dimensional  In¬ 
verse  Scattering  Methods,  M.  Cheney,  .1.11.  Hose,  and  H.  DcFacio 

We  relate  three  dimeii.sional  scattering  tlieoiy  for  the  time-independent  Schrodingi  r  equa¬ 
tion  without  spherical  symmetry  to  scattering  theory  for  the  plasma  wave  equation  (I’WLi). 
Wc  review  a  recent  inverse  scattering  method  for  the  I’WIO  and  find  the  corre.spond- 


ing  method  for  the  Sthiddinger  equation.  Wc  then  review  Newton’s  three  dimensional 
Marchenko  method  for  the  Schrddinger  ecpiation  and  transform  it  to  the  corresponding 


PWE  method.  The  resulting  time-domain  hyperbolic  method  clariries  the  role  of  causality 
in  Newton's  important  recent  work. 

A  Rigorous  Derivation  of  the  “Miracle^  Identity  of  Three  Dimensional  Inverse 
Scattering,  M.  Cheney 

The  large-energy  asymptotic  behavior  of  scattering  solutions  of  the  three  dimensional  time- 
dependent  Schrodingcr  equation  is  investigated.  The  second  term  of  the  expansion  leads 
to  the  “miracle"  of  Newton’s  three  (limcnsional  inverse  scattering  theory. 

Small  Gap  Flows,  E.O.  Tuck 

Tliis  report  consists  of  6vc  chapters,  namely: 

1.  Airfoils  in  Extreme  Ground  Effect,  pp.  2-6. 

2.  Automobile  Aerodynamics,  pp.  7-11. 

3.  Shallow-water  Waves  at  Di.scoiitimtities,  pp.  12-16. 

4.  Sliding  Sheets,  pp.  17-20. 

5.  Wiiig.s  Over  Water,  pp.  21-24. 

These  chapters  arc  independent  of  each  other,  although  linked  l)y  llic  “small-gap”  theme. 
C3ia))lcr8  1-4  were  originally  |>iesenled  in  the  forui  of  a  lecture  sei  it's  in  the  Department  (d 
Naval  Architecture  and  Olfshorc  Engineering,  University  of  California  at  lierkeley,  during 
February  and  March  1984.  Chapter  5  formed  the  basis  for  a  presentation  to  the  panel  115 
(Analytical  Ship- Wave  Relations]  of  the  Society  of  Naval  Architects  ;uid  Marine  Engineers, 
on  March  20th,  1984. 


Research  for  and  preparation  of  this  report  was  carried  out  during  a  period  spent  by 
the  author  in  the  Department  of  Mathematics,  Stanford  University,  January-May  1984. 


Discriminant,  Transmission  Coefficient,  and  stability  bands  of  Hill’s  equation, 
J.B.  Keller 

The  discriminant  A(fc^)  of  Hill’s  equation  is  shown  to  be  related  to  the  transmission  coef¬ 
ficient  of  one  period  of  the  potential  by  A(fc^)  [2/T'(fc)l  cos[fc7r  -|-0(A:)].  This 

result  is  used  to  find  the  boundaries  of  the  stability  bands. 

One  Hundred  Years  of  Diffraction  Theory,  J.B.  Keller 

The  development  of  diffraction  theory  in  the  last  100  years  is  discussed  from  a  personal 
viewpoint,  with  emphasis  on  the  geometrical  theory  of  diffraction.  First  some  c.arly  work 
of  Kirchhoff,  Rayleigh,  Somincifeld,  MacDonahl  and  others  is  mentioned  to  indicate  the 
state  of  the  field  in  the  1910’s.  Next  the  author’s  work  during  World  War  11  is  described. 
Then  the  considerations  that  led  hinr  to  the  geometrical  theory  of  diffraction  .are  expla’uied, 
and  the  defects  of  tliat  theory  are  outlined.  Fin.illy  the  advances  in  the  theory  since  its 
introduction,  which  have  remedied  many  of  these  rlefects,  arc  mentioned. 

Solition  Generation  and  Nonlinear  W\avc  Propagation,  J.B.  Keller 

Soliliou  gener.alion  by  v.iriou.s  means  is  described.  First  ex|»erimenl al  results  of  Wehausen 
ct  al  on  solilons  generated  by  a  ship  model  in  a  towing  lank  .ire  presenl«'d.  Then  T.Y.  Wu’s 
related  Boussinesq  system  of  equations  for  shallow  water  motion  with  a  moving  pressure 
distrubance  on  the  free  surface  is  introduced.  Numerical  solutions  of  this  system  by  Wu  and 
Wu  are  shown  to  compare  well  with  the  experimental  results.  Similar  numerical  results  on 


an  initial-boundary  value  problem  for  the  KdV  equation  by  Chu  et  al  are  presented  which 
also  yield  soliton  generation.  Then  Keener  and  Rinzcl’s  analysis  of  pulse  generation  in  the 
Fitzhugh-Naqumo  model  of  nerve  condutions  is  described.  Next  Whithain’s  modulation 
theory  of  nonlinear  wave  propagation  is  explained  and  the  problem  of  relating  its  results  to 
initial  and  boundary  data  is  mentioned.  Asymptotic  methods  for  solving  this  problem  for 
the  KdV  equation  are  described.  They  include  the  Lax-Levermore  theory  for  the  case  of 
small  dispersion,  its  extension  by  Venakides,  and  the  centered  simple  wave  solution  of  the 
modulation  equations  by  Gurevich  and  Pitaevskii.  Finally  the  theory  of  weakly  non-linear 
waves  of  Choquet-Bruhat  and  of  Hunter  and  Keller  is  described. 

Computers  and  Chaos  in  Mechanics,  J.B.  Keller 

Much  recent  progress  in  mechanics  is  related  to  the  use  of  computers  in  the  solution  of 
problems,  in  the  control  of  experiments,  and  in  the  an.alysis  of  data.  These  tools  have 
made  certain  kinds  of  calculations  and  mcasiircnicnts  easier.  But  they  have  also  revealed 
the  widespread  occurrence  of  chaotic  and  stocha.slic  behavior  of  mechanical  systems,  and 
have  shown  that  there  arc  certain  regularities  in  this  behavior.  The  understanding  of  the 
chaotic  behavior  and  its  rcgularilcs,  and  of  how  to  analyze  and  control  it,  arc  unsolved 
problems.  Many  theoretical  attacks  have  been  made  on  these  jirobicms,  but  they  have 
only  touched  the  surface. 

Sciniclassical  Mochaiiics,  J.B.  Keller 

Classical  mechanics  and  the  qutantum  conditions  of  Planck,  Bohr,  Sommcrfeld,  Wilson  iuid 
Einstein  arc  presented.  The  virtues  and  defects  of  this  "old  quantum  theory”  arc  pointed 


out.  Its  rcplnccnicnt  by  quantum  mechanics  is  described,  leading  to  the  Schrbclinger  equa¬ 
tion  for  the  wave  function  and  the  corresponding  energy  eigenvalues.  For  separable  sys¬ 
tems,  the  reduction  of  this  equation  to  ordinary  differential  equations  and  their  asymptotic 
solution  by  the  WKB  method  arc  described,  as  well  as  the  resulting  corrected  quantum 
conditions  with  integer  or  half-integer  quantum  numbers.  For  nonscparable  systems,  the 
analogous  asymptotic  solution  constructed  by  the  author  is  described,  together  with  the 
corrected  qiiantnni  conditions  to  which  it  leads.  Examples  of  the  use  of  these  conditions 


in  the  solution  of  eigenvalue  problems  arc  presented.  It  is  explained  that  difficulties  arise 
in  using  this  method  when  the  classical  motion  is  stochastic  or  chaotic.  Suggestions  for 
overcoming  these  difficulties  arc  mentioned. 
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